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A determination of the strange quark mass 



Roger Horsley 



1. The lattice approach 

Lattice methods allow, in principle, the complete 'ab initio' calculation of the fundamental 
parameters of QCD, such as quark masses. However quarks are not directly observable, being 
confined in hadrons and are thus not asymptotic states. So to determine their mass necessitates the 
use of a non-perturbative approach - such as lattice QCD. In this brief article, we report on our 
recent results for the strange quark mass for 2-flavour QCD in the MS-scheme at a scale of 2GeV, 
m" s (2GeV). Further details can be found in [p. 

1.1 Renormalisation group invariants 

Being confined, the mass of the quark, (M), needs to be defined by giving a scheme, 
and scale M, 

mf{M)=Z%{M)m™, (1.1) 

and thus we need to find both the bare quark mass and the renormalisation constant. An added 
complication is that the MS-scheme is a perturbative scheme, while more natural schemes which 
allow a non-perturbative definition of the renormalisation constants have to be used. It is thus 
convenient to first define a (non-unique) renormalisation group invariant (RGI) object, which is 
both scale and scheme independent by 

< = AZf {M)m y {M) = Z™m^ , (1.2) 

where, defining the coupling constant in the chosen scheme to be always g MS (i.e. expanding the 
fi y and Ym functions in terms of g" s ) we have 



[AZf (M)]- 1 = [2b g m (M) 2 ] ^ exp 




(1.3) 



The p ^ and Ym functions (with leading coefficients —bo, d,„o respectively) are known perturba- 
tively up to a certain order. In the MS scheme the first four coefficients are known, [§, ||], and this 
is also true for the RI'-MOM scheme [|], ||] (which is a suitable scheme for lattice applications). In 
Fig. [I] we show the results of solving eq. (L3) as a function of the scale M = jx and M = ji p for both 



the MS and RI'-MOM schemes respectively. We hope to use these (perturbative) results in a region 
where perturbation theory has converged. 2GeV corresponds to pL /A MS ~ 8, where it would appear 
that the expansion for the MS-scheme has converged; for the RI'-MOM scheme using a higher 
scale is safer (which is chosen in practice). However, when the RGI quantity has been determined 
we can then easily change from one scheme to another. Of course these scales are in units of A MS 
which is awkward to use: the standard 'unit' nowadays is the force scale ro. To convert to this unit, 
we use the result for roA Mi as given in pp. 

1.2 Chiral perturbation theory 

We have generated results for «/ = 2 degenerate sea quarks, together with a range of valence 
quark masses. Chiral perturbation theory, #PT, has been developed for this case, [[7], §]. We have 
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Figure 1: One-, two-, three- and four-loop results for [AZ% s (p)}- 1 and [hZ%- Mm '{lip)]' 1 in units of A* 



manipulated the structural form of this equation to give an ansatz of the form 

r mf = cf [(r m K+ ) 2 + {r m^f ~ {rm^f] 

+(cf-cf) [(rom K+ ) 2 + (r Q m K o) 2 ] {r m n+ f + \{cf + cf) [(r m K+ ) 2 + {rm&ff 

-{cf + c f){r m n+ f +cf( n) m n+ fHr m n+ ) 2 

-cf [{r Q m K+ ) 2 + {r m K o) 2 ] [{r m K+ ) 2 + {r m K o) 2 - (r m^+) 2 ] x 

In ({r Q m K+ ) 2 + (rom^o) 2 - {rotn K+ ) 2 ) +..., (1.4) 



and 



rom" 



f + cf{r m%) 2 +cf{r m ps ) 2 + cf {{r m%) 2 - 2{r m ps ) 2 ) ln(r m p ,) 2 . (1.5) 



{rom ps ) 



where m ps , m ps are the valence and sea pseudoscalar masses respectively (both using mass degen- 
erate quarks). The first term is the leading order, LO, result in ^PT while the remaining terms come 
from the next non-leading order, NLO, in #PT. We note that to NLO, we can determine cf and 
cf, i = b,c,d using mass degenerate quarks and then simply substitute them in eq. (Ob. 



1.3 The axial Ward identity 

Approaches to determining the quark mass on the lattice are to use the vector Ward identity, 
VWI (see e.g. |]9j]), where the bare quark mass is given in terms of the hopping parameter by 1 



m„ 



-\- 

2a \ K n 



1 



(1.6) 



or the axial Ward identity, AWI, which is the approach employed here. Imposing the AWI on the 
lattice for mass degenerate quarks, we have 



= 2m q ^> + 0{a 2 ) 



(1.7) 



This is valid for both valence and sea quarks. K^ c is defined for fixed /3 by the vanishing of the pseudoscalar mass, 

-qc ^ 



i.e. m ps {l^ c , k£.) = 0. k£, has been determined in 
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Figure 2: AZ£* (dashed line), Z A /Zf- M0M (empty squares) and Z*' (filled circles) for /3 = 5.40 (filled 
circles), together with a fit F{apL p ) = p\+ p2(an p ) 2 + p3/(afi p ) 2 . 



and 8$ and & are now the 0(a) improved 2 unrenormalised axial current and pseudoscalar density 
respectively and fn q is the AWI quark mass. So by forming two-point correlation functions with 9? 
in the usual way, this bare quark mass can be determined. We have found results for four /3-values: 
5.20, 5.25, 5.29, 5.40, each with three sea quark masses and a variety of valence quark masses, [jj|. 
Furthermore upon renormalisation we have that 

Zf (M) = f * . (1.8) 

As mentioned before, we use the Rf-MOM scheme, to determine Z A and z p ~ MCM in the chiral limit 3 . 
We now have all the components necessary to compute Z^' and hence rom". In Fig. || we show 
AZ%- U0M , Z A /Z^- M0M and their product, giving for J8 = 5.40. This should be independent of 
the scale (afl p ) at least for larger values. This seems to be the case, we make a phenomenological 
fit to account for residual {apLp) 2 effects. 

With Z", we can now find rora" and hence the ratio rora*'/ (rom ps ) 2 , using the values of ro/a 
given in [QJ. In Fig. || we plot this ratio (against (rom ps ) 2 ) for j8 = 5.40. Using eq. ( p~4] ) to eliminate 

in favour of ronif 3, /((romj f +) 2 + (rom^) 2 — (ro»%+) 2 ) in eq. ( |L~5| ) gives r^mf directly 4 to NLO 
in our fit function. 

We have restricted the quark masses to lie in the range (rom ps ) 2 < 5, which translates to 
m ps < 850 MeV, which is hopefully within the range of validity of low order #PT results. (Indeed 

2 The improvement term to the axial current, d p P together with improvement coefficient c A , [ |io| has been included. 
The mass improvement terms, together with their associated difference in improvement coefficients, b A , bp appear to be 
small and have been ignored here. 

3 For Z A (using the sea quarks only) we make a linear extrapolation in am q , while for Z p ~ M0M we must subtract 

out a pole in the quark mass, jllj], which occurs due to chiral symmetry breaking. We thus make a fit of the form 

{Zf- Mm )- l =A P +B P /am q . 
4t 



This is preferable to first determining c^ 1 and cj* 3 , i = 1,2,3 by using eq. (1.5 1 and then substituting in eq. fll.4|)) 
as the direct fit reduces the final error bar on rQmf J . 
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Figure 3: rom^ 1 /(rom ps ) 2 against (r(,m ps ) 2 , together with a fit using eq. (L5) for j3 = 5.40. Filled points 



represent valence quark results while unfilled points are the sea quark results. The dashed line (labelled 
'V^/Mjc') represents a fictitious particle composed of two strange quarks, which at LO ^PT is given from 



eq. (1.4) by yj (rom K +) 2 + (rom K o) 2 — (rom n +) 2 , while the dashed-dotted line (labelled 'm w ') representing a 
pion with mass degenerate u/d quark is given by r(,m K +. 



using ro/a, rather than their chirally extrapolated values for example, tends to give less variation in 
the ratio row"/ (r m ps ) 2 so we expect LO #PT to be markedly dominant.) Thus finally, for each 
j6 -value we have determined r^mf 1 and can now perform the last extrapolation to the continuum 
limit. 



2. Results 

Our derivation so far, although needing a secondary quantity such as ro/a for a unit, depends 
only on lattice quantities. Only at the last stage, with our direct fit did we need to give a physical 
scale to this unit. A popular choice is ro = 0.5 fm. However there are some uncertainties in this 
value; our derivation using the nucleon gave ro = 0.467 fm and so to give some idea of scale 
uncertainties, we shall consider both values. (The main change when changing the scale comes 



from the r s in eq. (1.4), as mf J <=c r , while changes in AZ™ S are only logarithmic.) 

Using the results from section [O] for [AZ^(2GeV)] _1 to convert mf to mf(2GeV) gives 
the results shown in Fig. ^ Also shown is an extrapolation to continuum limit. We finally find the 
result 

<(2GeV) = { 117 < 7 ! MeVft " r » = a5fm , (2.1, 
iV ' \ lll(6)MeV for r = 0.467fm ' 

where the error is statistical. This is to be compared to our previous result using the VWI, [g], 
which gave results of 126(5) MeV, 119(5) Me V for r$ = 0.5 fm and 0.467 fm respectively. We 
take a further systematic error on these results as being covered by the different ro values of about 
~ 6 MeV. Although the continuum extrapolation should be treated with caution, it does indicate 
that the strange quark mass for 2-flavour QCD lies in the region of 100 - 130MeV. 
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Figure 4: Results for mf i (2GeV) (filled circles) versus the chirally extrapolated values of (a/ro) 2 (as given 
in [^] ) together with a linear extrapolations to the continuum limit. For comparison, we also give our previous 
result using the VWI, ^ (open squares) and the ALPHA AWI determination from fll^], (open triangles). 
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